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Intermittency for branching walks with heavy tails. 

A. Getan * S. Molchanov J B. Vainberg 


Abstract 

Branching random walks on multidimensional lattice with heavy tails and a 
constant branching rate are considered. It is shown that under these conditions 
(heavy tails and constant rate), the front propagates exponentially fast, but the 
particles inside of the front are distributed very non-uniformly. The particles exhibit 
intermittent behavior in a large part of the region behind the front (i.e., the particles 
are concentrated only in very sparse spots there). The zone of non-intermittency 
(were particles are distributed relatively uniformly) extends with a power rate. This 
rate is found. 
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1 Introduction 

The mathematical stndy of branching processes goes back to the work of Galton and 
Watson [21] who were interested in the probabilities of long-term survival of family names. 
Later similar mathematical models were used to describe the evolution of a variety of 
biological populations, in genetics HQi 111121 da], and in the study of certain chemical and 
nuclear reactions nnun]. The branching processes (in particular, branching diffusions) 
play important role in the study of the evolution of various populations such as bacteria, 
cancer cells, carriers of a particular gene, etc., where each member of the population may 
die or produce offspring independently of the rest. 

In this paper, we describe the long-time behavior of a population in different regions 
of space when, in addition to branching, the members of the population move in space 
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(random migration). In fact, we will consider discrete problems on random walks on the 
lattice Z'^, not processes in bnt our results can be extended to the latter case (with 
some changes). 

The two main characteristics of the problem under consideration concern the rate of 
the branching and the probabilities of large jumps in Z'^. There is an extensive literature 
on the branching random walks with compactly supported or fast decaying branching 
rates, see [201 IH 13 l2l El El El El HEl |22l [18]. We assume that the branching rate v 
is constant (the total population grows exponentially with probability one in this case). 
Problems of this type (with constant rate v of splitting and zero mortality) originated 
from the KPP model, see the famous paper by Kolmogorov, Petrovski and Piskunov 
[15]. The spreading of a new advanced gene (in M*^) was studied in [15]. We consider 
the problem on the lattice (not in W^), but there is a much more essential difference. 
We allow the random walk to have long jumps. To be more exact, we consider random 
walks with heavy tails. This combination of heavy tails with the constant branching rate 
creates the following effect: the front of the population propagates exponentially fast, but 
the number of particles inside of the front is distributed very non-uniformly. The latter 
property is referred to as the intermittency for the number of particles, and is the main 
object of our study. 

Let us formulate the problem more precisely. Branching walks on the lattice Z'^ are 
considered. We assume that a particle located at ?/ G Z'^ at a given time t remains at 
the same point for a random exponentially distributed time (with parameter 1), and then 
jumps to a new position y + z with probability a{z). In addition, each particle located 
at ?/ G Z'^ splits with the rate udt into two particles located at the same point y. All the 
particles behave independently of each other and according to the same law. 

Obviously, 

a{z) >0, ^ a{z) = 1. (1) 

Z&d- 

We also assume that the distribution of the jumps is symmetric, i.e., a{z) = a{—z), which 
implies that L = L* in /^(Z'^). We assume that a has the following behavior at inhnity 

+ 1^1 ^ OO’ 

\Z \ \ ^\ 

where 

ao(i) > (5 > 0, 0 < a < 2. (3) 

In fact, later we will specify the asymptotics of a{z) in more detail (see fITOlB . 

Many essential characteristics of the random walk (or process in depend on the 
rate of decay of the distribution at inhnity. Conditions (|2]), ([3]) correspond to a walk with 
a heavy tail. The second moments do not exist in this case. The decay ([2]) with a > 2 
dehnes walks with a moderate tail. A process has a light tail if the distribution decays at 
inhnity so fast that its Fourier transform a{k) is analytic in k when |Im/c| < S with some 
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5 > 0. The second moments are well defined in the latter two cases. Thus the walk under 
consideration admits very long jumps whose probabilities are not very small. 

Let n{t, X, y) be the number of particles at a point x E a.t the moment t > 0 under 
condition that the process starts at initial moment t = 0 with a single particle located at 
a point y E i.e., n(0, x, y) = 5{x — y). 

Without duplication, the initial particle would perform the symmetric random walk 
X{t) with the generator 


{i^yf){y) = [f{y + ^) - f{y)]a{z). 

z£'L<i 


which acts on the space f (Z'^). 

If 1 / > 0 and irti = mi{t,x,y) is the expected value (the first moment) of the random 
variable n{t,x,y), then mi satisfies the relations: 


dmi 

dt 


(t, X, y) = Cymi{t, x,y) + v mi(f, x, y), 


t > 0. 


(4) 


mi{0,x,y) = S{x-y). 

In order to derive (jl]), we evaluate mi{t + At,x,y) by splitting the time interval 
(0, t + At) into two successive intervals of lengths At and t. Then 


mi{t + At, X, y) E a{z)Atmi {t, x,y + z) + 2vAtmi {t, x, y) (5) 

+ (1 — Af — vAt)mi{t, X, y). 

The first term on the right side of ([S]) is the sum of the probabilities to jump from y 
toy + z during the time At (these probabilities are a{z)At) multiplied by the expectation 
mi{t,x,y + z) for the number of particles at x when the walk starts at a single point 
y + z eZ/-. The second term describes the probability vAt of branching during the time 
At multiplied by the expected number of particles at x that are descendants of both the 
original and the new particles at y, which is 2mi{t, x, y). The last term is the contribution 
to the expectation of the number of particles at y from the event that the particle stays 
at X without branching and jumping during the time interval At. 

We subtract mi(t,x,y) from both sides of equation ([5]) above, divide by At, and pass 
to the limit At —)■ 0. This implies (jl]). 

Let us derive the equation for the second moment, m 2 (t,x,y) = E{in?{t,x,y)). We 
again consider the time interval (0, t + At) and split it in two successive intervals of lengths 
At and t. Then 

m 2 (t + At,x,y) E a{z)At m 2 (t, x,y + z) + V At E{ni + 722 )^ 

z&A 


+ {1 — At — 1 / At) m2(t, X, y). 


( 6 ) 


3 



Here the terms on the right of ([6]) are similar to the terms in ([5]), ni{t, x, y) is the number 
of particles in x that are descendants of the original particle, and n, 2 (t, x, y) is the number 
of particles in x that are descendants of the newly born particle. We use the fact that 
E{ni + 77 - 2 )^ = E{ni)‘^ + E(n 2 Y + 2E{ni)E{n2) = 2m2{t, x, y) + 2m\{t, x, y). Then we 
subtract m 2 (f, x, y) from both sides of (P), divide by At and pass to the limit as At —)■ 0. 
This implies: 


X, y) = {Cy + z/)m 2 (t, X, y) + 2v ml{t, x, y), (7) 

7772 ( 0 , X, 7 /) = 5{x - y). 

Denote by 7771 (t, x), 7772 (t, x) the hrst and second moments for the number of particles 
at the point x when the process starts at a single particle at the origin: 

777i(t, x) := 777i(t, X, 0), 7772(t, x) := 7772(t, X, 0) 

Obviously, functions mi(t, x, y), m 2 {t, x, y) depend on x — y, not on x and y separately, 

i.e., 

7771 (t,X, 7/) = 7771 (t,X - 2/,0), 7772(t,X, 7/) = 7772(t,X - 7/,0), 

and therefore, mi{t,x,y),m 2 it,x,y) are known as soon as 7?7i(t,x),7772(t,x) are found. 
Note also that the symmetry of a{z) implies that Cyf{x — y) = Ef{x — y), where 

Cf = C^f = [f{x + z)- f{x)]a{z) 

Hence mi{t,x), 7772 (t,x) satisfy the relations 

——(f, x) = (£ + 7/)777i(f, x), t > 0; 777i(0, x) = 5(x). (8) 

at 

^ (t, x) = (£ + Z/)7772(t, x) + 2z/ 777i(t, x), t > 0; 7772(0, x) = (5(x). (9) 

at 

The region in that separates the large and small values of mi{t,x) is called the 
front. To be more exact, we dehne the front E = E{t) as the boundary of the set 
{x : 777i(t, x) > 1}. The boundary consists of points x that have neighbors x + e-t, |e±| = 1 
or e± = 0, such that 777i(t, x + e+) > 1 , mi(t, x + e_) < 1 . 

The notion of intermittency (or intermittent random helds) is popular in natural sci¬ 
ences (astrophysics, biology, etc). From the qualitative point of view, intermittent random 
helds are distinguished by the formation of sparse spatial structures such as high peaks, 
clumps, patches, etc., giving the main contribution to the process in the medium. For 
instance, the magnetic held of the Sun is highly intermittent as almost all its energy is 
concentrated in the black spots, which cover only a very small part of the surface of the 
Sun. Many bio-populations also exhibit strong clumping (clustering). 
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Intermittency is a well developed non-uniformity. For physicists, the magnetic held of 
the Sun is intermittent since, say, 99% of its magnetic energy is concentrated on less than 
1% of the surface. For mathematicians, 0.1, 0.01 or 10“® are not necessarily small numbers, 
and a limiting process must be considered instead. The dehnition of intermittency based 
on the progressive growth of the statistical moments was proposed in the review [23], a 
more formal presentation can be found in |5]. In the the simplest form, a held n{t, x), x ^ 
is intermittent as f —)■ cxd on a non-decreasing family of sets D{t) if 

, En^(t,x) 
hm , ^ ^ = oo 

t^cx> (En{t, x)y 

uniformly in x E D{t). 

Let us illustrate this dehnition with the following example. Let n{t,x), x G be 
independent identically distributed r. v. and 

P{n(f,x) = 0 } = 1 - y, P{n{t,x) =t} = t>po. 

Then En{t,x) = po, i.e., the density of the population is constant in time, and 

Po = lim 
L —^oo 

However, 

En?{t, x) 

{En(t,x)r 

and similarly for , m > 2. Thus, the family of the helds n{t,x) is intermittent 

as t ^ oo. It is clear that the population n{t, x) is supported for large f on a subset of Z'^ 
with relative volume po/t —)■ 0 as f —)■ oo. The independence of n{t,x) is not important 
here and can be replaced by some kind of weak dependence (ergodicity), see [5|. 

The main two results of the present paper concern the propagation of the front of the 
branching random walk that starts at the origin, and the intermittency of the distribution 
of the particles on and behind the front. These results are proved under the following 
assumption that is a little more restrictive than (j^j). Namely, we assume that: 


(2L)-' 


(the law of the large numbers). 


d-\-e / . \ 




( 10 ) 


where i = zl\z\., 

ttj e ao(i) > 5 > 0, 

and e = lifQ; = l,e = 0 otherwise. 

Theorem 1.1. Let / fl^) hold. Then the following asymptotics holds for the points x G E{t) 
on the front E{t): 


x\ = [ao{x)t]<i+<^e'i+<^^{l + 0{t ^)), t ^ oo. 


( 11 ) 
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This theorem is an immediate consequence of the global limit theorem proved in [1] 
(see Theorem [2T] below and the remark after it). While the front propagates exponentially 
fast, the particles are distributed very non-uniformly on the front and at any exponential 
in time distance from the original particle. In fact, we will find the exact boundary for 
non-intermittency, and this boundary propagates with a power rate. Let 


2o -|- d 
a{a + d) 


( 12 ) 


The following statement will be proved in the next section. 

Theorem 1.2. Let / I7d|) hold. Then 

1) The ratio Is uniformly bounded in each ball |x| < BV when t —)■ oo, i.e., the 

random variable n is non-intermittent there. 

2) For each domain Tle{t) = {x : |a;| > e > t), we have —)■ oo uniformly 

in X & Tle{t) as t ^ oo, i.e., n is intermittent in fl^it). 


2 Proof of the main results. 

The function mi{t, x) with u = 0 will be denoted by p = p{t, x). Then 

^{t,x) = Cp{t,x) , t>0; p{0,x) = 5{x), (13) 

and 

mi{t,x) = e’'^p{t,x). (14) 

Consider the homogeneous (of order —d — a) distribution in that is equal to 

ao{x)\x\~^~^ when x ^ (compare with the first term of asymptotics ([2])). It was shown in 

[1] that the Fourier transform of this distribution is the homogeneous function —6o(d)|cr|“ 

in where a is the dual variable to a:, d = and 
’ ’dr 

QTT f 

6o(d) = —r(—a) cos — / ao(i)|(i, d)|"dS' 2 ; > 0. (15) 

2 

Here T is the gamma-function. 

The following global limit theorem for random walks with heavy tails obtained in [Ij 
is a key point in the proof of the results stated in the introduction. 

Theorem 2.1. Let ^ hold. Then 

(i) the following asymptotis holds for p 

xeZ'^, |x|+t^cx), (16) 
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were S{y) = > 0 is the stable density S = Sa,ao{y); which 

depends on a & ( 0 , 2 ) and Oq, and Bq is defined in 

U 1^1 ^ 1 ) then the previous statement can be specified as follows: 


p{t,x) 


ao{x) 

td/a 


(^)"“(l + o(l)) 

\x\ 


ao{x)t 

I ^ I (i-\-(y. 


(1 + 0 ( 1 )). 


(17) 


Remarks. 1. There is a misprint in the statement of the theorem in [T]: must be 

replaced by S'(^). 

2. This theorem and (ITTll imply Theorem 11.11 The remaining part of the section will 
be devoted to the proof of Theorem 11.21 

For a given function f = f{x), a: G denote by /(cr) the periodic in a function that 
is the Fourier series with coefficients f{x),x G i.e., 


m = Y 1 fio = 


(2ir) 




In particular, from fll3p it follows that 


p(t,x) = 


(2n) 


p(i, a)e‘^'’’‘)da = 


(2ir) 


^[a{a)-l]t+i{(T,x) 


a. 


(18) 


The following properties of a((T) follow immediately from properties of a{x)\ 


a(—cr) = a(a); — 1 < a{a) <1, 0 7 ^ a G T^. 


(19) 


The second part in (|T^ follows from ([T]) provided that for each a G T'^, cr 7 ^ 0, there is a 
point z E where 7 ^ 1 and a{z) 7 ^ 0. Such points z exist due to ffTOj) . 

We will need the following lemma. 

Lemma 2.2. Function p{t,x) is strictly positive for all x G f > 0. 

Proof. Denote by a„(a;) the convolution of n copies of a{x): 


an{x) := a(x) * a{x) * ... * a(a;), 


( 20 ) 


where a{x) * h{x) = — z) * b{z). 

We multiply both sides of fl20|) by and take the sum in a; G This implies 

On(cr) = [^(c)]*^- From here, ([T]), and flT^ it follows that |art(cr)| < 1 . This allows us to 
write p{t, x) as follows: 


p{t,x) 


1 

(2^ 


^la{a)-l]t+i{<x,x) _ 


[—7r,7r]^ 




nl 


Thus 
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Since a{x) > 0, all the convolutions a„(a;) are non-negative. Hence (j^ will imply the 
statement of the lemma if we show that a 2 {x) is strictly positive for all x G 
We have: 

02(3;) = a(x — z)a{z). ( 22 ) 

Here a > 0, and from (ITOll it follows that the terms in (|22ll are positive for each hxed 
a: if z is large enough. Thus 02 ( 0 ;) > 0, and the proof of the lemma is complete. 

□ 

Let us formulate a simplihed version of Theorem 12.11 that will be combined with fll4l) 
and will be easer to use than Theorem 12.11 

Definition. Functions a and b will be called equivalent, and it will be denoted by a x 6, 
if there exist two constants Ci and C2 such that, Cib < a < C26. 


Lemma 2.3. Let ^ hold. Then for arbitrary 02 > Oi > 0, the following relations hold 



(i) mi(f,x) X 1— when x > aif“. 

(23) 

(ii) 

mi(f, x) when x < 02 ^“, f > e > 0, 

(24) 


{Hi) mi(t,x) xl when x = 0, t < 1. 

(25) 

Remarks. 1) One can use any of the estimates fl25]l. (1241) in the intermediate zone 
Oit^ < X < 02 ^“, t > e > 0. The right-hand sides of these estimates are equivalent in 
this intermediate zone. 

2) From the last two relations it follows that the following estimate holds for all t > 0: 


mi(f,x) X (t-|-1) when x < 02 ^“. 

(26) 


Proof. Due to flTT)) . it is enough to prove the lemma when i/ = 0 and mi is replaced by 

p. 

The last statement of the lemma follows immediately from fflSj) . In order to prove the 
hrst two statements, we split the region x G > 0, into three subregions, Ui, U 2 , and 
U 3 , where Ui is the region dehned by the inequality |a;| > Af^ with A so large that the 
remainder term in ffT7|l is less than 1/2. Note that |a;| > Afz implies that x ^ 0. Thus 
|a:| > 1 in Hi since x is a point on the lattice. Hence flT7)l implies (12^ with oi = A. Let 
U 2 be dehned by the inequalities |x| < AT^, |x| -|- t > H, where B is chosen so large that 
the remainder term in (1T6|1 is less than 1/2 when \x\+t > B. Since -^ < H is bounded in 

t OL —— 

U 2 and function S{y) is positive and continuous, it follows that ^(-x) in flTB]) has upper 

t 

and lower positive bounds in f/ 2 . Thus flTB]) implies that p{t,x) x in H 2 . 








Now consider the region 

t /3 = {|x| < Af^, |x| + t < -B}\{(t,x) : X = 0 , t < 1 }. 

Region U 3 is bonnded, and t > 6 > 0 there. Since p is continnons, from Lemma 12.21 it 
follows that p has positive lower and npper bonnds on t/ 3 . Similar bonnds are valid for 
the fnnction t~'^. Thns p{t,x) x t~a on 1/2 U ^ 3 ; the second statement of Lemma 
12.31 is proved with 02 = A. Estimates (12^ . were proved with oi = 02 = A. Their 
validity with arbitrary oi, 02 follows from the eqnivalency relation stated in the remark 
above. 

□ 

Proof of Theorem 11.21 Since rrii is the Green fnnction for the operator C + u, the 
Dnhamel principle implies that the solntion m 2 of fl27D has the form 

m 2 (t, x) = mi(f, x) + 2i/ / mi{t — s,x — z)m‘l{s, z)ds. (27) 

We will start the proof of the theorem with the second statement (abont the intermit- 
tency). Since mi,m 2 > 0, we have 


m 2 (t, x) > 2i/ 



mi{t — s, X — z)ml{s, z)ds, 


zSZ"*: h|<5 


1 

ta 


(28) 


where x G Ge(t), t ^ 00 . We apply estimate fl23l) to the hrst factor nnder the snmmation 
sign and the estimate fIMl) to the second factor. Taking into acconnt that t — s < 1 and 
|x — ^1 X |x| in fl28l) . we obtain that 


m2(t,x) >C ^ 


it-i ^ 1 


\X 


— s “ e^’^^ds > C 

\d+a — 


f*t — 

t°‘ 2ut 7 

s “ e ds. 


n-i 


X 


\d+a 


Hence 


2 i't 


t-|x|‘^+" 

On the other hand, from fl23D it follows that 

Ct^ 


/ ^ C'e ^ / X 

m 2 {t,x) > — -, X G t —>■ 00 . 


mf{t,x) < . '6 ^^ , X G Oe(t), t ^ 00 . 


X 


\2{d+a) 


Hence, if x G t ^ 00 , then 


^ CO a. 7 


ml{t,x) I'i+d/a ^ 


—)■ 00 . 
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The second statement of the theorem is proved. 

Let us prove the hrst statement. Note that the ball |x| < BP is located behind the 
front, and Lemma f2.3l implies that mi —)■ cxo in this ball as t —)■ cxo. Hence fl27j) implies 
that the hrst statement of the theorem will be proved (and the proof of the theorem will 
be complete) as soon as we show that 

1 V 

I :=—57-r / 7 mi(t — s, X — z)ml(s, z) ds < C, \x\ < BP, t — >■ oo. (29) 

In order to estimate the left hand-side of fl2^ . we split in (12^ into four sets 
separated by the two spheres: \z\ = and \x — z\ = {t — s)^. Let Pi = Pi(s),P 2 = 
P 2 {x,t — s) be (bounded) sets of points z & located inside or at the boundary of the 
hrst (respectively, second) sphere dehned above, i.e.. 

Pi = {z e Z"^ : |z| < S“}, P 2 = {z e : \x — z\ < {t — s)^} 

Denote by Dij = Dij{s, t, x) the following sets of points z eZ^ with 0 < s < t, x E Z^: 

Dll is the the set of points z E Z'^ located inside or at the boundary of both spheres, 
i.e., Pii = Pi Pi P 2 . 

D 22 is the the set of points z E Z'^ located outside of both spheres, i.e., P 22 = 
Z''\(PlU^2). 

P 12 is the set of points z E Z'^ located outside of the hrst sphere, but inside of the 
second one or on its boundary, i.e., D 12 = (Z'^ \ Pi) IJP 2 . 

P 21 is the the set of points z eZ'^ located inside of the hrst sphere or on its boundary, 
but outside of the second one, i.e., P 21 = (^'^ \ P 2 ) U -Pi- 

Respectively, / can be written as / = where lij is dehned by fl2^ with the 

summation extended oved Dij instead of Z'^. We are going to estimate each of the terms 
lij when PI < BP, t —)■ 00 . 

1) Estimate on I 12 . We will obtain this estimate separately for 2 ta < |x| < BP and 
for |x| < 2 ta. Note that the two balls Pi and P 2 are separated in the hrst case (when 
0 < s < t) and they may intersect each other in the second case. Consider the hrst case 
of PI > 2P. Then \mi{t,x)\ x (due to (l23lL and relations fl26l) . (I23|) hold for 

the hrst and second factors under the summation sign in fl2^ . respectively. Thus 


h 2 < 


C\xf^ 




2Mt+s)cis 


s e 


[(t-s) + l]^P|2C+“) ’ 


2t°‘ < PI < BP, t —)■ cx). (30) 


Since p — z| < (t — s)“ < in P 12 , inequality p| > 2 t^ implies that p — ^| < ip|, 
and therefore p| > |p| in fl5U]l . Hence we can replace zhy x there. After that, summation 
in fl30|) is applied to z-independent terms. Hence, the summation sign can be replaced 
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by a factor k that estimates the number of terms in the sum from above. Obviously, 
K < C{A‘^ + 1), where A is the radius of the ball P 2 , i.e., n < C[(t — s)'^ + 1]) ^ind 


I (i-s + i)l 


< 


[ s^e'^^ds < C 2 < 00 

7o 


(31) 


when 2 ta < |x| < BP, t ^ 00 . 

Now let |x| < 2 ta. Then \ml{t,x)\ x for f > 1 (due to (Ell)), and (13(1 must be 

t a 

replaced by 


/12 < 


Ct 


^ rt 


o 2 ut 


2Gi?12 O 


s^e‘ 


ds 


s + l)^\z\^A+A 


\x\ <2t°‘, t ^ 00 . 


(32) 


We split the right-hand side above into two parts I' -|- J" by writing the interval (0, t) 
as the union of (0,f/2) and {t/2,t). Since D 12 does not contain points of the hrst ball Pi, 
2 ; 7 ^: 0 there. We have \z\ > 1 for all other points of the lattice. Thus we can replace \z\ 
in r by 1 and replace summation by factor k. This leads to the estimate 


/' < 


C-tw p/2 [(t - s)l -f /•*/" 


— ^ 2 i't 




< 
— 




ds, 


where the right hand side decays exponentially as t ^ 00 . Thus I' < C < 00 . 

Now consider J". In this case, t/2 < s < t implies that \z\ > > (t/2)“, and we 

can replace \z\ by (t/2)“ in fl32l) . Then we replace the summation in z by the factor k as 
before. This leads to the following estimate: 


I" < 


CP 


^2ut-i^2{d+a) I a 


h/2 


< C < (50, 


t —)■ CX). 


Together with the boundedness of /' and (l30|) . this proves that 


1\2 < C < 00 , |x| < BP, t ^ 00 . 


2) Estimate on In. The two balls Pi, P 2 do not intersect each other when |a:| > 
2P (and 0 < s < t). Hence we may assume that |x| < 2 ta. Then fl24p implies that 
\m‘l{t,x)\ X t~~, t —)■ cx). We apply fEB]) to the factors under the summation sign in 
(E^ and obtain that 


/n < 


Ct 


^ rt 




o 2 ut 


2gr>ii ((: - s -M)«(s -M) 


t —CX). 
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We replace here the summation sign by the factor k introduced above and use the 
estimate 

K, 


{t - s + 1 )'< 




(t — s + 1 ) ' 


This leads to 


/ii< 




d'^ds 


Jo (s + 1 )^ 


< C < CX), t —)• CX). 


( 33 ) 


3 ) Estimate on l22- Estimate fl 2 ^ can be applied to the factors under the summation 
sign in fl 2 ^ . Thus 


I 22 < 




l{t,x) Jo ^ \x — 


( 34 ) 


Consider hrst the case when |x| < \t^^°‘- Then the inequalities 
|;2| > S“, |x — z| > (t — s)“ , 0 < s < t, X, 


( 35 ) 


imply that \z\ >/ 3 > 0 iff = l. Using the homogeneity arguments, one can easily obtain 
that inequalities fl 3 ^ with an arbitrary t > 0 imply that 


i.e., the latter estimate holds in ^22- Now we can replace \z\ in fl 34 D by and use 

relation flMll for mi. This leads to 


-^22 < 


(j^2d/c 


^2ut-j-2(d+a) / a 


E 

ZGD 22 


{t — s)s'^e''^^~^^^ds 


X 


I d-\-a 


|x| < 

' ' 2 


t —)■ 00. 


Note that |x - 2;| > 1 when 2; G -D22 and that Y.Z&D22 |x-^V+" ^ ^z(izd\{x} The 

latter series converges and does not depend on x. Hence 


/■* 1 
I22 < I ~ s)s‘^e’^^^~^^'^ds < C < 00, |x| < t —)• 00. ( 36 ) 


'0 


Now let us estimate I22 when < |x| < BE. We split the region D22 into two, 

namely, £>22^ = D22 n {z : [z] > -y) and = D22 Cl {z : |z| < -y}- Then I22 < h + h, 
where Ji and I2 are the right-hand side in fl 3 T)) with D22 replaced by ^ D^22 respectively. 

In order to evaluate Ji, we replace \z\ in flMjl by |x|/ 2 , use convergence of the series 
YI1Z&D22 jW2p+“’ estimate fl 25 ]l for mi. This leads to 


C\x\^(d+a) 
— i2g2ti/ 



{t - s)s^ 

^ 12(d+Q) 


e‘'^*+^^ds 


Ct-^e-^^ 



s)e'^^ds < C < 00 
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when < |x| < ^ t —f oo. 

Let us estimate J 2 . The inequality |^| < -^ implies that \x — z\ > We use this in 
(l3T|) together with estimate fl2^ for mi. Then we have 


l2< 


C\x\(^+») 


{t 


s)s^e'^H 




zGD. 


( 2 ) 


j|2(d+a) 


)ds, 


1 .l/c 

2 


< |a:| < BB, t —)■ cx). 


Note that \z\ > Sa in D 22 and \z\ > 1 ii z 0. Hence the sum in the formula above does 
not exceed C(1 + s)~^^. We also can use there that ^ (due to |a;| < BB). 

Hence 

^j-dja. pt -1 

I 2 <—:— / (1 + s) ^(t — ds, < \x\ < BB, f —)■ 00 . 

e*'" Jo 2 

One can easily check that the right-hand side above is bounded. This together with the 
boundedness of Ji and (1361) implies that 


I 22 B C < 00 , |x| < BB, t ^ 00 . 


4) Estimate on Bi- From Lemma 12.31 it follows that 


Bi < 


C 


m\{t, x) 


Z&D 21 


{t-s) 


\x — z 


d-\-OL 


{s + iy 




First, let us estimate Bi when |x| < 2 ta. Then mi can be estimated using fl2TD . There 
is also x-independent constant O' < 00 such that ^ ^ details in the 

subsection on B 2 )- Thus 


Bi < 


CB- 




o2vt 


I / N -iCl 

B (s-I-1) 




CB> 


oVt 


r (t-s) 

lo (s + 1)^ 


■e^®ds < C < 00 


when |x| < 2 B, t —>■ 00 . 

Consider now the case of 2fx < |x| < BB. Then estimate (1231) can be applied to mi, 
and therefore 


L 21 < 


C\xyid+a) 

^2g2i/t 



{t — ^ 

OJ OjO . 

X - zy+°‘{s + 1)1^ 


(37) 


Since \z\ < < t^ in D 21 and we assume that |x| > 2 ta, we have |x —> |x|/2 in the 

integrand above. Thus we can replace \x — z\ by |x|/2 in (157|) . After that, the summation 
sign can be replaced by the number ki of terms in the sum. Obviously, ki < C{Af + 1), 
where Ai is the radius of the hrst ball Pi, i.e., Ki < C[sa -|- 1], and 


B^ < 


C|x 


d+a 


^2gi/t 


{t — s)e 

(s + l) 


US 

x^ds < - 


{t — s)e^ 




(s + l) 


ds < C < 00 


d+1 — 
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when < |a;| < t —)■ cxd. Thus I 21 is bounded when |a;| < BV , t ^ 00. Together 

with the boundedness of all other lij, this completes the proof of the theorem. 


□ 
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